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We study the distribution of eigenvalues for the Green operator occurring in the scattering of electromag- 
netic waves by an arbitrarily shaped dielectric medium. It is revealed that the totality of eigenvalues 
(counting multiplicities) can be enumerated as a sequence {A s } _,,N < Ko, with only two possible accu- 
mulation points {0,-1/2}, and the following spectral series converges: Z _i |A S | |1 + 2A S | 4 < +oo. 



1. INTRODUCTION 

In [Ref. 1, Chap. 4], we studied the analytic properties of the Born equation for electromagnetic scattering: 

E(r')e- iklr - r ' 1 



fE(r')e- lK,r - r ' , , 
4n|r-r'| dr = E ™ {r) > reV - (1) 



Here, the bounded open set VgK 3 represents the dielectric volume occupied by a homogeneous dielectric medium with 
complex susceptibility % - e r - 1, where e r is the relative dielectric permittivity, n - s/e^ the complex refractive index. 
We require that the boundary surface dV be smooth and the exterior volume R 3 \ (V U dV) be connected, but the shape 
of V can be otherwise arbitrary. The (real part of the) complex-valued incident field E{ nc (r)e lckt ,r e V represents a light 
beam oscillating in the fixed angular frequency cj-ck, where c is the speed of light. Here, the time-harmonic factor 
e lckt will be conveniently suppressed hereafter. The dielectric response (total electric field) inside the volume V is given 
by E(r),r e V. We abbreviate Eq.[l]as (I - x c S)E - £ mc , with / being the identity operator, and w the Green operator. 

It is a physical requirement that both E and E mc be square integrable and divergence-free. Therefore, we will pose 
the Born equation (7 - x<g)E = E inc (Eq.Q) on a physical Hilbert space <t>(V;C 3 ) = Cl(C°°(V;C 3 )nker(V-)nL 2 (V;C 3 )), 
which is the L 2 -closure of the totality of smooth, divergence-free and square-integrable complex- valued vector fields, 
[cf. Ref.Q, Chap. 4] The physical spectrum rj*(<#) of the Green operator <g : <P(V;C 3 ) — <D(V;C 3 ) is the totality of com- 
plex numbers A e C that forbid a bounded inverse of the operator XI — *& : <1>(V;C 3 ) — ► <1>(V;C 3 ). Practically speak- 
ing, a point in the physical spectrum points to an "optical resonance mode" where the energy amplification ratio 
fff v \E(r)\ 2 d 3 r/fff v \E inc (r)\ 2 d 3 r goes without bound. 

We have proved, in [Ref. 1, Chap. 4], the following theorem regarding the topological structure of the physical spec- 
trum ct*C#). 



Theorem 1.1 (Compact Polynomial and Optical Resonance) The non-compact Green operator <£ : <1>(V;C 3 ) — ► 
%7 i <J>(V;C ) is polynomially compact, with minimal polynomial 'Sil + 2 ( £)I2. 

Moreover, cr®(w) is the union of the continuous spectrum afCS) - {0,-1/2} with the point spectrum o-pCiS) that con- 
tains countably many eigenvalues. Each eigenvalue has a strictly negative imaginary part, 1 and is associated with a 
finite-dimensional eigenspace. The continuous spectrum crfCS) - {0,-1/2} forms the only possible accumulation points of 
eigenvalues. 

The shape-independent singularity -1/2 e afCS) corresponds to a universal optical resonance at susceptibility % - -2, 
i.e. relative permittivity e r --l. | 

In this work, we will provide more detailed information regarding the spectral structure of <£ : <1>(V;C 3 ) — ► <I>(V;C 3 ) 
by proving the following theorem. 



'Electronic address: yajunz@math.princeton.edu 
The connectedness of B 3 \(Vu dV) is conducive to ruling out real-valued eigenvalues, as well as confirming that e a *("$). Even when R \(Vu dV ) 
is not connected, the smooth dielectric boundary dV is strong enough to ensure the compactness of c S(i + 2 e £):<&(V;C ) — » S>(V;C ). 



Theorem 1.2 (Hilbert-Schmidt Polynomial and Spectral Series) The operator <S{I + 2 ( £) 2 : d>(V;C 3 ) — d>(V;C 3 ) is 
of Hilbert- Schmidt type, and we have a convergent spectral series 2 

£ |A| 2 |l + 2A| 4 <+oo, 

where the sum respects multiplicities in eigenvalues. | 

Theorem 11.21 refines Theorem 1 1.11 by providing more quantitative details about the structure of the physical spectrum 
cr*(w). While 'Sit + 2^) is already compact, it takes another factor of (7 + 2^) to yield a Hilbert-Schmidt operator 
<£(I + 2 C S) 2 . While both the point {0} and {-1/2} can function as accumulation points in the spectrum a^CS), the two 
distinct powers |A| 2 and |1 + 2A| 4 occurring in the convergent spectral series indicate different levels of eigenvalue 
aggregation surrounding the two accumulation points. 

This article is organized as follows: §2 sketches the proof of polynomial compactness for 'S : <1>(V;C 3 ) — ► <J>(V;C 3 ), 
thereby recapitulating the key ingredients of Theorem 11.11 and providing some analytic motivations for Theorem ll.21 §3 
gives a detailed proof of Theorem [1.21 by extending the thoughts debriefed in §2; §4 provides some concrete examples 
that verify Theorem [L2] and discusses its physical interpretation as well as possible applications. 

2. POLYNOMIAL COMPACTNESS OF THE GREEN OPERATOR <S 

The Born equation (Eq. Q} is a strongly singular integral equation (with integral kernel behavior 0(\r-r'\~ 3 ) and 
the Green operator ^ is not compact. Nonetheless, the Green operator ^ can be naturally decomposed into a sum of a 
Hermitian operator and a compact operator as 'S - CS - f) + f. Here, the Hermitian operator 'S - f is given by 

E{r ' ] -dV 



Iv 4:ii\r — r'\ 



E(r')(e- iklr - r ' l -l) 3 



(d-f)E)(.r):=V V- fff 

and the compact operator f is defined as 

(fE)(r) := V x V x fff r^TL d V. (2) 

JJJv 4n\r-r'\ 

Physically speaking, the operator <S — f represents the long-wavelength limit (2n/k — * +00) of electromagnetic scat- 
tering, which is just electrostatic induction; the operator 7 represents "dynamic corrections" on top of the electro- 
static picture. Mathematically speaking, the "static induction" operator ^ - 7 is Hermitian, satisfying the inequal- 
ity > (E,(<S -f)E)v a -(E,E)v [Ref. 1, p. 188]; the "dynamic correction" operator 7 is weakly singular, with only 
OQr-rT 1 ) divergence in the integral kernel [Ref. 1, p. 176], hence is qualified as a Hilbert-Schmidt operator, which is 
a special type of compact operator. 

As the difference between the quadratic polynomial <£ + 2 ( S 2 and the expression ('S - 7) + 2(^1 - f) 2 is evidently a 
compact operator, we only need to verify the following proposition before claiming the compactness of if + 2 ( S 2 . 

Proposition 2.1 The expression (<# - 7) + 2(<£ - f) 2 : <1>(V ;C 3 ) — ► ®(V ;C 3 ) defines a compact operator. 
Proof (Sketch) A key observation is that, the "static induction" <S - 7 maps any E e ®(V ;C 3 ) to 



m-f)E)(r): 



JJd 



n'-E(r') 
'dv 4ji\r-r'\ 



as' 



which is the gradient of a harmonic function, reminiscent of an electrostatic field. The bulk behavior (CS-f)E)(r),r e V 
("electrostatic field") is fully determined by the boundary behavior of its normal component n ■({ ( S- f)E)(r),r e dV 
("boundary surface charge"), and this one-to-one correspondence is also robust in both directions [Ref. 2, p. 252]. Tech- 
nically speaking, according to the boundary trace theorem and the robustness of Neumann boundary problems, there 
are two finite positive constants C\ and C2 such that 

Ci||(# - f)E\\ L 2 (v . c3) < ||»- ((<# - y)E)\\ H -i,2 (av . c) < C 2 ||(» " f)E\\ L 2 (V . C 3 } . 



The convergence remains true even when the exterior volume R \ (V u dV) is not connected. 



The compactness of the polynomial {0 — j)+ 2{ c £ - f ) 2 on <J>(V ;C 3 ) is then evident from the boundary integral repre- 
sentation 



n ■((<$- f)E)(r) + 2n ■ {{<$ - f) 2 E)(r) = -2 E n' «<S - f)E)(r')(n ■ V) * - dS', re dV 

UdV 4n\r-r'\ 



(3) 



where the integral kernel (»• V)(47r|r-r'|) _1 = Odr-r'l -1 ) is weakly singular (see [Ref. 3, p. 48] or [Ref 4]), hence 
induces a compact operator (cf. [Ref. 4] or [Ref. 5, p. 124]) on H~ 1/2 (dV;C). Here, the choice of the coefficient 2 in the 
quadratic term of { c S—f) + 2(<S - f ) 2 is critical, in that 



lim(n-V)<£f n'-((<3-i)E)(r')— -dS' 

£-0+ JJdv 4ji\r-en-r'\ 

E n' ■ {{<£ - f)E)(r')(n • V) - 1 dS' + Ji» • ((» - f )£)( 
JJdv 4n\r-r'\ 2 



r), redV. 

(Here, the limit is interpreted in the sense of boundary trace.) In other words, the construction of the quadratic polyno- 
mial ensures that the non-compact ingredients cancel out. | 

In the next short lemma, we will derive a volume integral analog of Eq.[3j after introducing the notations for electro- 
static Green functions. The Green function Go(r,r') with Dirichlet boundary condition is given by the unique solution 
to 

V 2 Gz>(r,r') = -<5(r-r'), r,r' eV; G D (r,r') = 0, reV,r'edV. 

The Green function thus defined automatically honors reciprocal symmetry Go(r,r') - Go(r',r) [Ref. 6, p. 40]. Custom- 
arily, the Green function GN(r,r') with Neumann boundary condition is prescribed as a solution to 

V 2 Gjv(r,r') = -<5(r-r'), r,r' eV; (n' ■V')G N (r,r')= - - \ - - , reV,r'edV. 

Hereafter, we will impose the constraint of reciprocal symmetry to the Green function GN(r,r') = GN(r',r) (see [Ref. 01 
or [Ref. 6, p. 40] for such feasibility). 



Lemma 2.1 Define 



Q(r,r'):=G D (r,r') + G N (r,r')- - * r,r'eV, 

2n\r-r'\ 



which is a harmonic function with respect to both r and r', then we have the volume integral representation 

((<£-f)E)(r) + 2((<£-f) 2 E)(r) = JJJ VV' (r,r')-((^-f)£)(r')dV, reV, (4) 

where VV'g(r,r') is a 3 x 3 matrix filled with mixed second-order derivatives ofg(r,r'). 
Proof By properties of the Green functions, we have 

[[[ VV' (r,r')-((^-f)£X''')dV = Vc£f fl (r,r')i»'-((»-f)E)(r')dS' 
JJJV JJdV 

= VE G N (ryW.^-f)E)(rW-2vE »'<f-W\ S , 
JJdv JJdv 4n\r-r'\ 

= (d - f)E)(r') + 2{{<£ - f) 2 E)(r), reV, 

as claimed in Eq.[4j | 

From the above arguments, it is straightforward to check that the integral kernel K\(r,r') :- VV'g(r,r') is Hermitian 
K^(r,r') - K\{r' ,r), where an asterisk denotes the conjugate transpose of a 3x3 matrix. What we have proved in 
Proposition 12. 1 1 is that the integral kernel VV'g(r,r') induces a compact linear operator, unlike the strongly singular 
integral kernel associated with the non-compact operator C S. In the next section, we will show that the following integral 
kernel (where products represent matrix multiplications) 

fff VV" fl (r,r")V"V' (r",r')dV' = JJJ K^ry'W^r" ,r')d 3 r" 
is square-integrable, so it induces a Hilbert-Schmidt operator, just as the "dynamic correction" operator f . 



3. HILBERT-SCHMIDT POLYNOMIAL IN THE GREEN OPERATOR & 

According to the Schur-Weyl inequality (see [Refs.0,131, also [Ref. 10, p. 8]), we have 



, X IA| 2 H + 2A[ 4 <||W + 2'#) 2 ||2:=WEll^ + 2^) 2 e s ||2 2(yc 3 ) 

\ A£(7*(W) V 8=1 

where {e s \s - 1,2,...} is any complete set of orthonormal basis for the Hilbert space <1>(V;C 3 ), so the Hilbert-Schmidt 
bound ||^(/ + 2^) 2 1| 2 < +oo would entail the convergence of the spectral series in question. Using the inequalities 
II AS || 2 < l|A|| 2 ||B|| and ||BA|| 2 < l|A|| 2 ||S|| [Ref. H p. 218], together with the facts that \YS-j || < 1, ||/ + <#-y|| < 1 and 
II f II - II f II 2 < +°°, we may deduce the inequality 

||(/ + 2<#) 2 c#||2 

a\2/ 



||(7 + 2^-2f) z (^-f)||2 



„"Vv2„-, 



ft\„",2 



+ ||f + 4[(^-f)f + f(^-f)](/ + ^-f) + 4(^-yrf + 4(7 + ^- 7)7^ + 47(^-7)7 + 4f z (^-fr + 4f d || 2 



.<Mi2x 



<||(7 + 2^-27) z (^-f)||2 + |l7ll2(13+12||f||+4||f|| z ). 



(5) 



Therefore, the major task in Theorem ll.2l boils down to an analysis of the action of (7 + 2^ - 2f) 2 on the range of (<S - 7), 
which is represented by the integral kernel 



fff W"g(r,r")V'V fl (r",r')dV'. 



Before establishing the Hilbert-Schmidt bound for the operator (7 + 2^ - 2f) 2 { ( £ - 7) in Proposition 13. II we explore the 
boundary behavior of the integral kernel g(r,r') in the lemma below. 

Lemma 3.1 7<br any f £ 77" (dV;C), we have the following limits in the sense of boundary trace: 



(«'-V) 



1 1 

■+ ■ 



2n\r-r'\ § dV dS 



1 1 

(n-V) — : - + - 



2n\r-r'\ § dV dS\ 



AS'; 
AS', 



lim <ff f(r')(n'V%(r-en,r f )dS'=-E f(r') 

c^O+JJdV JJdV 

lim(»V)<££ g(r-en,r')f(r')AS'=-E f(r') 

£— 0+ JJdV JJdV 

where r e dV. 

Proof We first prove (*) by computing 

lim <fcb f(r')(n'-V')Q(r-en,r')AS' 

f-0+ JJdV 

= limc£f f(r')(n'-V')G D (r-en,r')AS'- ■= * ff f(r')AS'-2 lim & f (r')(»' • V)— -AS' 

f-o+JJay § dV ASJJdv n-o+JJdv 4ji\r-en-r\ 

= -f(r f )- -= * ff f Xr') AS'- 2 lim E f(r')(n'-V)— -AS' 

§ dv ASJJdv E~o+JJdV 4n\r-En-r'\ 



(*) 

(* *) 



-«ff f(r')\(n 
JJdV 



■V) .* ..+■ ' 



2n\r-r'\ § dV AS 



AS', reV, 



where we have invoked the solution to the Dirichlet boundary value problem using Go , and the Neumann boundary 
condition for Gn- To tackle (* *), we perform the following analysis: 



lim(»V)<ff Q(r-en,r')f{r')AS' 

£-0+ JJdV 

- lim(»-V)(Q) 

£-0+ JJd 



G N (r-en,r')- 



2n\r - en - r' \ 

; av f(r")AS"i 



f(r')AS' 



= lim(n-V)<££ G N (r-en,r') f(r) 

£-0+ JJdV [ 

-2 lim E f(r')(n-V)— -AS' 

£-o+JJav 4n\r - en - r' \ 

= -$ fir') 
JJdV 



ffdV dS 



dS'+ mv ' jo lim(n-V)<H> G N (r-en,r')AS' 

J 3fdv AS £-0+ JJdV 



(n-V) 



1 1 

■ + ■ 



2n\r-r'\ § dv AS 



AS', redV, 



s 

= -1 + 



where we have relied on the homogeneous Dirichlet boundary condition for Go , the solution to the Neumann boundary 
value problem in terms of Gn, and the fact that 

lim(nV)<jf£ G N (r-£n,r')dS' 

£-0+ JJav 

lim lim(re-V)$ G N (r - en,r' -8n')dS'+ lim lim(re-V)$ [G N (r- en,r')-G N (r-En,r' -8n')]dS' 
>-o+£-o+ JJdv 5-0+ £-0+ JJgv 

lim lim(re-V)c/£ [G N (r-£n,r')-G N (r-En,r'-8n')]dS' = 0. 

>-0+ £-0+ JJgv 

Here, the last line can be justified as follows. For sufficiently small 8 > 0, we may define the boundary layer of thickness 
8 as Vr s) := {r e V|dist(r,3V) < 8], so that the function dist(r,d\0 := min r / ea y i — r'\ is smooth in r e V~ udV,~ , satisfying 
(re- V)dist(r,dV) = -1. Thus, choosing v' as the outward normal of the smooth domain VZy we have 

lim lim(re-V)c££ [G N (r-en,r')-G N (r-en,r' -8n')]dS' 
5-0+ £-o+ JJav 

= lim lim(re-V)c/£ [dist(r' ,dV)(v' -V')G N (r- En,r')-G N (r- £n,r')(v' -V')dist(r' ,dV)]dS' 

5-0+ £-0+ JJdV,- 

(8) 

= lim lim (re-V)fff [dist(r',dV)V 2 G N (r-en,r')-G N (r-£n,r')\/' 2 dist(r',dV)]d 3 r' 

5-0+ £-0+ JJJv,- 

(o) 

= - lim lim(re-V)dist(r-£re,dV) = l. I 

5-0+ £-0+ 

Now, for fixed r' £ dV, the boundary trace (re ■ V)g(r,r'), as a function of r e 5V, has merely a weak singular- 
ity Odr-r'l -1 ), hence (re-V)g(r,r') £ #"~ 1/2 (dV;C). As a result, in the independent variable r, the harmonic vector 
field Vg(r,r') e L 2 (V;C 3 ) is square integrable, i.e. g(r,r') e W r ' (V;C) for all boundary points r' e dV. Furthermore, 
su.p r t £dv fffv |Vg(r,r')l 2 d 3 r is finite. The boundary trace mapping from the Sobolev space W 1,2 (V;C) toH 1/2 (dV;C) then 
leads to Q(r,r') e HJ. /2 (dV;C),r' e dV. Moreover, the expression fff v \Vg(r,r')\ 2 d 3 r defines a subharmonic function in r', 
which may only attain its maximum at the boundary r' edV. Therefore, the condition sup r / e5 y fff v \V$(r,r')\ 2 d a r < +00 
entails the square integrability of V g(r,r') e L 2 (V ;C 3 ),r' eV. As a result, we always have Q(r,r')eH]! 2 (dV;C), no mat- 
ter the point r' is at the boundary r' e dV or in the interior r' e V . 

Proposition 3.1 The operator (I + 2<i- 2f) 2 { ( S - f ) : ®(V; C 3 ) — ► d>( V; C 3 ) is of Hilbert-Schmidt type. 

Proof Without loss of generality, we may pick the complete orthonormal basis set {e s |s= 1,2,...} c<J>(V;C 3 ) so that one 
of its subset {f s \s - 1,2,...} exhausts all the eigenvectors subordinate to the non-zero eigenvalues of the polynomially 
compact Hermitian operator & -f : <!>( V;C 3 ) — ► <t>(V;C 3 ). The bound on the operator norm \\& - f|| < 1 then naturally 
leads to 



||(/ + 2^-2f) 2 (^-f)|| 2 <W£l|(/ + 2^-2f) 2 /- s ||2 2(vc3) . 

V s = l 

In the formula above, each f s is the gradient of a harmonic function, so the following integral representations hold for 
FeClspan{/- s |s=l,2,...}: 

((/ + 2^-2y)F)(r)= fff Wg(r,r')F(r')d 3 r' = fff £i(r,r')F(r')dV, 
((7 + 2^-2f) 2 .F)(r)= fff VV'<5(r,r')F(r')d 3 r' = HI K 2 (r,r')F(r')d 3 r' , 
where 

6(r,r'):= fff V'g(r,r")-V" S (r",r')d 3 r" 

evidently satisfies the harmonic equations V 2 &(r,r') - and V' 2 &(r,r') - for r,r' e V. 

Now, extending the action of the integral kernel K2(r,r') on [f s \s = 1,2,. . .} c ®(V ;C 3 ) to a complete orthonormal basis 
set of L 2 (V ;C 3 ), and using the Parseval identity on L 2 (V ;C 3 ), we can show that 

2 |ia + 2^-2f) 2 /- s ||^ 2(V;C3) = J^|j^Tr[^|(^,r-')^ 2 (A-,r-^]d 3 ^|d 3 ^, (6) 



where "Tr" denotes the trace of a 3x3 matrix. To justify the equality in the above formula, we note that 
fff v K2(r,r')F(r')d 3 r' represents the gradient of a harmonic function for whatever square-integrable input F e 
L 2 (V;C 3 ), so every eigenvector F A e L 2 (V;C 3 ) satisfying fff v K 2 (r,r')Fx(r')d 3 r' = AF A (r), A ? must belong to the func- 
tion space <I>(V;C 3 ). In other words, the extension of the orthonormal basis set leaves the Hilbert-Schmidt norm intact. 
We may go on to cast the volume integral representation of &(r,r') into a surface integral (more precisely, a canonical 
pairing [Ref. 2, p. 206] between g(r,r") e H^(dV;C) and (»" • V")g(r",r') eH:,j /2 (dV;C)) in two ways 



<6(rJ)=& 0(r,r")(»"-V")fl(r",r')dS"= ff (r",r / )(»"- V")fl(r,r")dS", 

JJdV JJdV 



r,r'eV. 



Clearly, the reciprocal symmetry g(r\,r2) - Q(r 2 , ri) entails the result &(r,r') - ®(r',r). Interpreting the expression 
&(r,r'),r e V,r' e dV as a boundary trace (denoted by the limit notation "lim E ^ +" as before), we may use Lemma lcTTt *) 
to deduce 



<8(r,r')= lim ff g(r,r")(n" -V")g(r",r' -en')dS" 

c-0+ JJdV 



ff Q(r,r") 
JJdV 



(n -V ) 



2n\r'-r"\ § gv dS' \ 



dS", reV,r'edV. 



Then, employing the harmonic equations for <3(r,r'), we may convert the double volume integral on the right-hand side 
of Eq.[6]to a double surface integral in the following fashion: 



!!! Iff! Tr[£ 2 *(r,r')£ 2 (r,r')]dVld 3 r= [[[ £ J [[[ £ [(».V)(p-V')<S(r,r')] 2 dV I d 3 r 
JJJvUJJv J JJJv u e{ ex , ey , ez }[JjJV vs{ex!eyie:s] J 

= ((( E Iff [(»-V)(S(r,r')][(«-V)(n'-V')0(r,r')]dS'ld ; 

= «£f Iff [(»-V)<5(r,r')][(n'-V')0(r,r')]dS'ldS. 

J Jsv I J Jay ) 



(7) 

Here, the integrands in the last line are understood in terms of boundary trace, expressible as a specific case of 
Lemma l3. ![ **): 

(n-V)<S(r,r') = lim( M -V)©(r-en,r') 



JJd\ 



(»-V) 



1 1 

■ + ■ 



Wl 2jr|r-r"| $ 5y dS \ 



(n"-V") 



2;r|r'-r''| # flV dS'J 



dS", r,r'e<3V. 



Judging from the surface integral above, the singular behavior of («• V)©(r,r') is comparable to the convolution of two 
OOr-r'l -1 ) integral kernels on the boundary surface, which results in the short distance asymptotics (n-V^ir,^) - 
0(\n\r-r'\). Likewise, by reciprocal symmetry ®(r,r')= <3(r',r), we have 



(n'-v')<3(r,r') = (n'-v')<S(r» = ff 

JJdV 



(ra'-V) 



1 1 

■ + ■ 



2n\r'-r"\ § av dS> 



(n"-V") 



1 1 

■ + ■ 



2n\r-r"\ § dv dS 



dS", 



which is again a surface integral kernel of order 0(ln|r-r'|). As the logarithmic singularity is square integrable, the 
surface integral 



ff [(i»-V)0(r,r')][(j»'-V')0(r,r')]dS' 

JJdV 



is finite for every r e dV, it is then evident that the double surface integral in Eq. [7] converges. Hence, we have estab- 
lished the Hilbert-Schmidt bound ||(/ + 2<#-2y) 2 (#-y)|| 2 < +oo. | 

Remark In the course of this proof, the expression in Eq. [7] provides a (theoretically) computable upper bound for 
\\(I + 2<S - 2f) 2 (<S - f)\\ 2 . For an arbitrary shape VdR 3 , this eventually boils down to a quadruple surface integral of 
purely geometric quantities: 

V)<S(r,r')][(n'-V')<3(r,r')]dS'U<S 



JJdV {JJdV 

= ff dSxff dS 2 ff dS 3 ff dS 4 Y\ 

JJdV JJdV JJdV JJdV j-\ 



n j-( r j - r j+lmod4) 1 

2n \ r j ~ r j+l mod 4l 3 §dV AS J 



which could be practically challenging in numerical evaluations. In the current work, we will not further discuss the 
practical computation of the surface integral in Eq. [7J except for a specific example in 34-11 concerning a spherical 
boundary surface dV. 

It might appear as an interesting fact that the degree of the Hilbert-Schmidt polynomial deg( ( S(I + 2 ( S) 2 ) = 3 = dimlR 3 
equals the dimension of the space where the dielectric volume V resides. This is not a pure coincidence, which we will 
explain later in 94.31 | 

4. EXAMPLES AND DISCUSSIONS 

4.1. Hilbert-Schmidt Qualifications in Mie Scattering 

In the Mie scattering scenario, the dielectric medium occupies a spherical volume V - 0(0, R) with radius R. Using 
spherical harmonics Y( m (8,cf>), we may verify Theorem [L2] with brute force, and compare the geometric representation 
of the Hilbert-Schmidt norm (Eq.Q with an algebraic approach. 

Proposition 4.1 The operator ( S{I + 2 ( S) 2 : d>(O(0,i?);C 3 ) — ► <P(O(0,i?);C 3 ) is of Hilbert-Schmidt type, while the operator 
C S(I + 2 ( S) : ®(O(0,R);C 3 ) — ► <J>(O(0,.R);C 3 ) is compact, without being a Hilbert-Schmidt operator. 
Moreover, we have the identity 

f \\(I + 2y-W 2 fs\\ 2 moi0R) . c3) = $ \& [(n-V)6(r,r')][(»'-V')«(r > r')]dS'ldS = f — i-j 

s =l L (U(U,K),L ) JJ d O(0,R) (JJdO(0,R) J fT' 1 (2£+l) d 

for any complete orthonormal basis {f s \s - 1,2,...} of Q>(O(0,R);C 3 ). 

Proof It would suffice to check the first two claims on the operator polynomials (I+2 ( S-2j) 2 { ( S-f) and (I +2^ -2f)( ( S - 
y). According to [RefQ, P- 190], we have the spectral decomposition: 

(<$-f)E = -t t '^"'fW'" , V£e«>(O(0,E);C 3 ), 

e=i m =-e 2£+l 

with the complete orthonormal basis set {// ro (r) = (£R 2e+1 Y ll2 V(\r\ e Y em {8,(p))\£ = l,2,...;m= [-£,£]nZ] for Cl((<#- 
f)$(O(0,E);C 3 )). 

By definition of the Hilbert-Schmidt norm, we may compute 



/ °° £ 2 
\\(i + 2<S-2r) 2 (<£-i)\h = J I ^rj^g < +oo, 

so (/+ 2<S - 2f ) 2 (<£ - f) : <P(O(0,i?);C 3 ) — <P(O(0,i?);C 3 ) is a Hilbert-Schmidt operator. 
Meanwhile, the spectral decomposition 

(I + 2<S-2f)CS-f)E = -f: £ £{fe ™; E \ V I em ' VEe<P(O(0,R);C 3 ) 

tells us that the operator (/ + 2<# - 2f)(<sl - f) : <D(O(0,i?);C 3 ) — ► d>(O(0,E);C 3 ) is the uniform limit of a sequence of 
finite-rank operators, hence compact. However, this compact operator fails the Hilbert-Schmidt criterion because 



||(7 + 2^-2f)(^-f)ll 2 = jE 7^7-^3 = +cx >" 



tl (2£ + 1) 3 



This stark contrast reveals that the additional factor (/ + 2^) plays an indispensable role in turning a compact operator 
W + 2^) into a Hilbert-Schmidt operator ■#(/ + 2 C S) 2 . 
From an algebraic perspective, we doubtlessly have 

OO OO (! OO 1 

E ik/ + 2^ - ar) 2 /^ ii£ a(0(0 ^ >;C3) = E E iitf+2^-2f) 2 /*ji£ 2 c8) = £ j^-^- 

s=l (=\m=-l e=1 (4Z + 1) 
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However, we may also compute the Hilbert-Schmidt norm by explicitly evaluating the geometric integral (Eq. 0, with 
the help of spherical harmonics. Concretely, we have 



4n\r-r'\ 



-2RnV- 



2R 



n-(r-r') 1 



4n\r-r'\ 4n\r-r'\ 



'|3 






e=0m 



according to the simple algebra 2Rn ■(r-r') = r-r-2r-r' + r'-r' = \r - r'\ 2 for r,r' e dV - dO(0,R). This leads to the 
calculation 



(wV)©(r 



4n z Ud 



dO(0,R) 
p2 1^ 1^ 



n-(r-r") 1 



|r-r"| 3 2R 2 
1 



n"-(r"-r) 1 
|r"-r'| 3 ~23 



dS" 



and consequently Eq.[7]can be computed as 



Y;(6,<l>)Y,(e',<l>') = in' • V')(8(r,r'), 



2 em y 



ff Iff [(n-V)<8(ry)][(n'-V')<8(r,r')]aS'\dS=T—± 

JJdO(0,R) (JJdO(0,R) I e=1 (2£ + 



7C(3) 



+ iy 8 



1. 



While it becomes impractical to algebraically enumerate all the eigenvectors of the operator 'S - f for non-spherical 
shapes, we may still fall back on the geometric integral representation (Eq.[7]l to deduce a finite Hilbert-Schmidt norm 
from the bounded curvatures of the smooth boundary dV. | 

Remark For the sphere with radius R , the Green function with Dirichlet boundary condition can be developed into the 
following spherical harmonic series [Ref 6, p. 65]: 



G D (r,r') = 



1 . + £ £ Y Lw>V Y ?J e '>^ 



min(|r|,|r'|) 



|r|VTmax(|rUr'|) 



R 



ie+\ 



4;rmax(|r|,|r'|) 4nR £ lm f_ { (2£ + l)max(|r|,|r'|) [ [max(|r|,|r'|) J 
1 1 °° < \r\ e \r'\ e Y; m {8,<p)Y em W',<p') 

~ 4n\r-r'\ 4nR £i„,t/ {2£+l)R 2e+1 

Meanwhile, the symmetrized Green function with Neumann boundary condition has an explicit series representa- 
tion [Ref.S p. 144]: 



G N <r,r') = 



1 f f VJ^Y^Jff^ 

47rmax(|r|,|r'|) f-' lm tL e (2£ + l)max(|r|,|r'|) 



min(|r|,|r'|) 



max(|r|,|r'|) 



+ l|r|Vrmax(|r|,|r'|) 
~£ R 2t,+1 



- t e + 1 \r\ e \r'\ e Y; m {e,<p)Y p W',<p') 



^ + E E 



4n\r-r'\ e=lm= 
It is then clear that 

Q(r,r'):=GD(r,r') + G N (r,r')- 



(2£+l)R 2e+1 



1 1_~ ' \rf\r'\%* m W,<p)Y em W',<P') 

2n\r-r'\~ 4nR + ' fciJh-t 



£{2£+l)R 2e+1 
provides yet another means to solve the eigenvalue problem 

((I + 2<S-2f)F x )(r) = [ff VV'Q(r,r')-F A (r')d 3 r' = XF A (r), F A e<P(O(0,i?);C 3 ). 

JJJO(0,R) 

Furthermore, the series representation of g(r,r') can be summed into 



1 rl OO I 

0(r,r')+-i-= £ E 

4nR Jo e = lm tL 



loo / f-WlrrYTJOrfW.JO'rf) 



(2£+l)R 2 ^ 



r t\r\\r'\ r' 



Jo l||r| R 



2 i_/ 



Ir'l 



d£ 
47ri?i 



-In- 



1 -LTL + \JL 

1 R 2 + \\r\ R 2 |r'| 



Thus, there is a mild logarithmic divergence g(r,r') - 0(ln|r - r'\) if at least one of r and r' is situated on the spherical 
boundary dO(0,R). Overall, the singular part of VV'0(r,r') has order 0(\r-r'\~ 2 ). Such asymptotic behavior is not 
a trend specific to spherical geometry, but a generic property descending from the curvature bounds of the boundary 
surface dV (see g4.3t . 

Owing to the bound ||^-y|| < 1/2 in the spherical case, the Hilbert-Schmidt norm estimate in Eq. [5] does allow a 
modest improvement in the case of Mie scattering: 

||(/ + 2<#) 2 <#|| 2 



< 



ii(/+2^-2fm-f)ii 2 

+ ||f + 4[(<#-f)f + f(^-f)](/ + ^-f) + 4(^-f) 2 f + 4(/ + ^-f)f 2 + 4f(^-7)f + 4f 2 C#-f) 2 + 4f 3 ||2 



OO £2 

I 7^7-rT5 + Hfll2(6 + 7||f||+4||f|| 2 ). 

Here, the infinite series in question evaluates to a small number 

g <^ = 703) 3105) _^ =(00g21+)2 
^i(2/+l) 5 32 128 192 

and the upper bounds of H7II2 and ||y|| will be given in the next subsection, in a more general geometric context. | 



4.2. Bounds of Spectral Series in Arbitrarily Shaped Dielectrics 

In this subsection, we will describe the bound estimates of ||f|| and ||f H2 for dielectric media with arbitrary geometry. 
Combining these estimates with Eq. [7j we will arrive at a theoretical upper bound for the convergent spectral series 

Xa e ^)W 2 H + 2A| 4 . 

In [Ref. 1, pp. 261-272], we established the following bound for ||f||. 

Proposition 4.2 For a bounded open dielectric volume V with smooth boundary dV, we have 

>2/3 



max |ReA| < 1+ ||Ref|| < 1 + mhJ — Ik 3 fff d 3 r) ,3\/2kR v l, 



A£(7 Q (») 



4:71 JJJv ' 5n { JJJv 
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and 



- min ImA= max |ImA| < ||Imf || < min-{ — /// d"V, — [h* III dV ,-kR v 
Ae<r*(W) Aecr*(#) [ 4:71 JJJv 5n \ JJJv I 8 

where Ry :- min reK 3 macx. r i £ vudv \ r ' ~ r \ is the minimal radius of all the circumscribed spheres. | 

Now, we will turn to provide an upper bound for the Hilbert-Schmidt norm ||y||2- 
Proposition 4.3 For arbitrarily shaped bounded open dielectric volume V, we have 

,2/S 3 

where Ry is the minimal radius of all the circumscribed spheres. 

To prove Proposition [4J3J we need some preparations. First, we turn Eq.|2]into the following integral representation 

(fE)(r):=fff f(r,r')£(r')dV 

where the 3x3 matrix T(r,r ! ) is given by 



ll;> < mm<{ ^ \k 3 fff d 3 r] ,^(kR v f 



f(r,r') = 



4n\r-r'\ 



-ik\r-r'\ 



l-e- ik ^- r ^{l + ik\r-r'\) 



k 2 \r-r'\ 2 



+ 



^ r - r \l + l k\r-r'\)-l _ lW _,, 

3 ^^7^2 < e -i) 



rr 
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Here, 

^^T _ r ~ r I r-r j _ , A .T\T _ A ^T - <*T 



, =(rr 1 ) 1 =rr'rr 1 
\r-r'\[\r-r'\) 

is a projection operator. Then, to facilitate further analysis, we may use the breakdown f-f^ + f^ where 

(y„.E)(r) := JJJ f„(r,r')£(r')d 3 r' = jfjfj^ ^~^r [e- ikir -'h + (1 - e-'^^' 1 )^] £(r')dV 
and 

Lemma 4.1 T/ie following results hold: 

3k 4 
Tr[ff(ry)f < (ry)]= ltoa|r _ ry|a , 

T rf * r /v f , m 3|l- e -' fe ' r - r ''(l+^|r-r f |)| 2 ^ 3£ 4 

Proof We note that for a,b e C, there is an algebraic identity 

Tr[(a1l + 6rr T )*(al + 6rr T )] = Tr[|a| 2 l + 2Re(a*6)rr T + |6| 2 rr T ] = 3|a| 2 + 2Re(a*6)+|6| 2 . 

In particular, setting a - e ~ lklr ~ r ' } and b - l-e"^ |r_r| , we have 3|a| 2 + 2Re(a*6)+ \b\ 2 = 3; setting a = 1 and 6 = -3, 
we have 3\a\ 2 +2Re(a*b)+\b\ 2 = 6. For x > 0, we have |e" -l-ix| < x 2 /2, which confirms that Tr[f *(r,r')t u (r,r')] < 

Tr[fj(r,r')f,(r,r')]/2. | 

Now we may turn to consider the functional 



w W//Il//I^F dV K 



so that we have the Hilbert-Schmidt bounds 

V3k 2 



£{£ Tr[f » (r ' r ' )f » (r ' r,)]d3r J d3r ^ /m ' 



and 

V6k 2 



£{///v Tr[f « (r ' r ' )f « (r ' r,)]dV } d8 ^-S- /[71 - 



As a result of triangle inequality for the Hilbert-Schmidt norm, we obtain 

2 + || f Qt | || 2 < IiVl 

Oil 



Lemma 4.2 We have 

(rrr ^ 
^£ d3r 

and consequently, we have \\fW2 £ (2/jt) (& 3 jQfjfy d 3 r) 
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Proof Using the Hardy-Littlewood-Sobolev inequality rRefs. fl34l3l with Lieb's sharp constant (see [Ref. lal . as well as 
[Ref. 16, Theorem 4.3]): 



If ( 

\Ju d Jk 



f(r'Mr) Ad d 

■d r d r 



«d \r-r'\' 



<7T 



Yf d—m \ 
ml2 ly 2 ' 

rw-f) 



r(f) i 
rw) 



II f \\ L 2d/(2d-m) (Jjd .Q II h || jidl&d-m) ( K d .Q , 



where d = 3,m = 2,— = -; fir) = h(r) = \ ' 

2d-m 2 \0, r$V 



4/3 



we may arrive at the following estimate 

Accordingly, the estimates on ||y|| 2 derives from 1.997+ = (2^+ v / 6)(v / 27r) 2/3 /8 < 2. | 

Remark From this lemma, we have I[O(0,R)]/R 2 < 7.022+, which may be compared to the exact value I[O(0,R)]/R 2 - 
2n obtained in the next lemma. | 



Lemma 4.3 We have 



I[V]- 



\M v e^dW 

r 3 47r|g| 



= v//l 



d 3 g 



and in particular, we have the identity I[O(0,R)] - 2nR 2 for all R > 0, and \\fW2 £ ^(kRv) 2 - 
Proof We start with the Fourier inversion formula 



fff 1 dV- 1 fff 27l2 Wv el9 ' r ' d3r ' c -i«r di 

n)v\r-r'\ 2 (2n? JJJ& \q\ 



where the convolution kernel satisfies 



lim [[[ 



r 

lim / 

-0+J0 



+oo e -eW sin (| g || r |) 2 7T 2 



inn li! -^-e i " r d 3 r = 4 7 rlim / ,""7"' " d|r| = 

|r| 2 £ -0+Jo \q\\r\ \q\ 



Then, applying the Parseval-Plancherel identity, we obtain the Fourier representation of I[V] as claimed. 
For V = O(0,R), we have 



1 



iq . r , , 3 , sin(|g|B)-|g|Bcos(|g|fl) A „ z M\q\R) 



V=O(0,R) 



\q\ d 



■ AnR A 



\q\R 



where j\{x) is the spherical Bessel function of first order. Consequently, we may derive I[O(0,R)] - 2nR 2 from the 
identity J +oo j 2 (x)/xdx = 1/4, which is a specific case of the Weber-Schafheitlin integral [Ref.[l7|, PP- 402-404]. We may 
check this against the spherical harmonic expansion 

~ / Y; m W,<l>)Y em W',cp')\rmn(\r\,\r'\)^ 
4n L L 



\r-r f \ 2 



r 0m tl/(2^ + l)max(|r|,|r'|) 



max(|r|,|r'|) 



R 1 r R 



which integrates to 

//I ( o,«) {//l ( o,«) ^P d V ' } d3 ^ (47r)2 £ c^TT) X" IX 

00 

= (4n) 2 R 4 Y, 



[min(|r|,|r'|)] 2 ' , , |2j| ., ., 



t 2(2£+l)(2£ + 3) 



[max(|r|,|r'|)] 



= (2nR 2 ) 2 . 



MY12/+2 



|r'| z d|r'|V|r| z d|r 



Suppose that we have chosen the origin r - of the coordinate system to coincide with the center of the circumscribed 
sphere of radius Ry, then we have the relations V c O(0, Ry), I\V] < I[O(0,Rv)], and this eventually establishes 

||f||2<(2v / 3 + v / 6K^v) 2 /4<|(^y) 2 . | 
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So far, we have verified Proposition 14.31 in its entirety. We may check the reasonability of the estimates in Propo- 
sitions 14.21 and 14.31 by revisiting to the spectral series for Mie scattering on a dielectric with radius R in the short- 
wavelength limit kR —>■ +oo: 

£ |A| 2 |l + 2A| 4 <||7l||||2f|| 4 = 0(* 8 fl 8 ). 

Here, the 0(k 8 R 8 ) behavior can also be recovered from the left-hand side. In the short-wavelength limit, the dominant 
contributions to the spectral series £ Aeo .< S( c| ) |A| 2 |l + 2A| 4 would arise from the eigenvalues associated with Y( m modes 
for £ < kR. Each such eigenvalue has order O(kR) in magnitude, and there are about 0(k 2 R 2 ) terms of them, counting 
multiplicities. Therefore, the spectral series ZA e(7 *(c|)|A| 2 |l + 2A| 4 = 0(k 8 R 8 ) is roughly the sum of 0(k 2 R 2 ) addends, 
each with magnitude 0(k 6 R G ). 

Even for non-spherical dielectrics, the consistency of Z^ E£r O(c|)IA| 2 |l + 2A| 4 < ||f || 2 ||2y|| 4 in the short-wavelength limit 
is ready to check. As H^-fH < 1 remains bounded irrespective of wavelength, we may pretend that the spectra o r ("$) 
and cr®(f) are close to each other in the short-wavelength limit, so that we have the following intuitive estimates: 



E IA| 2 |1 + 2A| 4 : 
Aeo-®(#) 



£ |A'| 2 |2A'| 4 <l|2f|| 4 £ |A'| 



'|2 



A'ecr*(f) 



A'eo-*(f) 



l|2f|| 4 . 



4.3. Comparison with Acoustic Scattering and Low-Dimensional Electromagnetic Scattering 

The analytic peculiarity of electromagnetic scattering problems can be better appreciated if we contrast it to the 
scattering problems concerning scalar acoustic waves in spatial dimensions d - 1,2,3. 
For example, the acoustic scattering equation in three-dimensional space 



((I-X^)u)(r):=u(r)-x 



E 



u(r')k 2 e- iklr - r ' 1 

4n\r-r'\ 



d r =u inc (r), 



eL 2 (V;C) 



involves a perturbation of the identity operator 7 by a Hilbert-Schmidt operator -%<€ :L 2 (V;C) — ► L 2 (V;C), and thus 
has significantly simplified spectral properties as compared to the Born equation for light scattering. The integral 



equation (7-;^)u = Wj nc is a Fredholm equation with convergent spectral series £ 



Ae(r(^) 



|A| Z = I 



AetjpCif) 



|Ap < +oo. Such 



an integral equation admits relatively simple expansion of (/-^Sf) -1 for any 11% t crC^) [see Ref. 10, Chap. 5]. 

For the scattering of electromagnetic waves, we lack compactness in the Green operator ^ : $(V;C 3 ) — ► <1>(V;C 3 ), and 
it takes a quadratic polynomial ( £{I + 2 ( S) :<1>(V;C 3 ) — ► <1>(V;C 3 ) to recover compactness, a cubic polynomial ( S(I + 2 ( S) 2 : 
<t>(V;C 3 ) — ► <I>(V;C 3 ) to qualify as a Hilbert-Schmidt operator. 

Heuristically, we can interpret the compactness of (<£ -f) + 2{<£ - f) 2 : <t>(V ;C 3 ) — ► <1>(V ;C 3 ) in terms of an integral 
kernel VV'g(r,r') with 0(\r-r'\~ 2 ) asymptotic behavior, which mitigates the strong singularity of order 0(\r-r'\~ 3 ) 
associated with the non-compact operator ( S. Here, the harmonic function Q{r,r') develops singularity as the points r,r' 
approach the dielectric boundary dV, and their distance approaches zero \r-r'\ — ► + , so we may sketch the 0(\r-r'\~ 2 ) 
behavior of the integral kernel VV'g(r,r') by focusing on the near-boundary cases. 

We choose r' e dV as a boundary point and use T r i(dV) to denote the tangent plane of dV at the point r' e dV, and 
write n' for the outward unit normal vector at r' . We situate the interior point r-r'- en' e V so close to the boundary 
point r' that the distance e :- \r-r'\ is negligible as compared to both the characteristic linear dimension £y of the 
volume V and the principal radii of curvature at r' e dV. The local "surface charge density" (n* • V»)g(r*,r') for a 



Integral Equation 


Spectral Series 


Operator 


Integral 
Kernel 


Compact 


Hilbert- 
Schmidt 


Electromagnetic Scattering 

a- X <£)E = E ine 


<^:<D(V;C 3 ) — <D(V;C 3 ) 

^Ae ( T*(«) |;l|2|1 + 2;l|4<+00 


» 


0(|r-r'r 3 ) 


No 


No 


^(7 + 2^) 


0(|r-r'r 2 ) 


Yes 


No 


»(7 + 2#) 2 


OCIr-r'r 1 ) 


Yes 


Yes 


Acoustic Scattering 

(7 - t&)u = u inc 


^L^C) — L 2 (V*;C) 


"% 


OCIr-r'r 1 ) 


Yes 


Yes 



Table I A comparison of various operators arising from three-dimensional electromagnetic and acoustic scattering problems, in terms 
of their asymptotic behavior in the integral kernel, their compactness, and their Hilbert-Schmidt boundedness. 
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Integral Equation 


Spectral Series 


Operator 


Compact 


Hilbert- 
Schmidt 


Electromagnetic Scattering 

(i- x <£(d) )E=E . nc 


I |A| 2 |l + 2A| 2W - 1) <+oo 


cg(d) 


1 


1 


# d) (/ + 2# d) ) 


1,2,3 


1,2 


# d) (/ + 2# rf) ) 2 


1,2,3 


1,2,3 


Acoustic Scattering 


I |A| 2 <+oo 


<£«) 


1,2,3 


1,2,3 



Table II A comparison for the dimension-dependence of the spectral structure in electromagnetic and acoustic scattering problems. 
The last two columns are filled with spatial dimensions d where the operator in question is compact, or Hilbert-Schmidt. 

boundary point r» £ dV near r' has a singular part ~ xy/|r, - r'\ where K r i is a quantity comparable to the curvature of 
the surface at point r' . The "electrostatic potential" g(r,r') thus scales asymptotically as 



Q(r,r 



■>-£. 



K r /d 2 (r» -r') 



W)nO(r>,e v ) 4n\r t -r'\\/E 2 + \r*-r'\ 



^2 Jo 



2 v / i 2 + p 2 2 A 



K r /dp 



K r / , 2^y 

— In . 

2 \r-r'\ 



This rough estimate hints at the 0(\r- r'\~ 2 ) bound of the integral kernel VV'g(r,r'), which hearkens back to the 
spherical example computed in the Remark to Proposition [47T] Intuitively, one may say that the 0(|r-r'|~ 2 ) bound for 
the volume integral kernel descends from the boundary integral kernel with Odr-r'l -1 ) behavior in Proposition ^. II 
Such an intuition can be made rigorous for arbitrary boundary geometry, once we work out the compactness proof of 
Proposition 12. 1 1 in fuller detail in the volume integral setting [Ref. 1, pp. 203-214, pp. 232-246], where the Calderon- 
Zygmund cancellation conditions [Ref. 18, p. 306 and p. 324] explain the mollification of the strong singularity 0(|r-r|" 3 ) 
in the operator polynomial (& — f) + 2(<S - y) 2 . 

Consequently, in view of the integral kernel VV'g(r,r') with 0(\r- r'\~ 2 ) behavior, along with the identity [Ref. [l^, 
p. 117] 



JJJs 



r3 2tt 2 \ 



I 2 2tt 2 | 



_/|2 



dV' : 



4n\r-r'\ ' 



we may envision that the integral kernel corresponding to (I + 2<£ - 2f) 2 (<S - f) would have Odr-r'l -1 ) asymptotic 
behavior, which is square integrable. This gives an intuitive understanding of Theorem 1 1.2 1 

Next, we show that the analytic structure of wave scattering not only is affected by degrees of freedom (vector wave 
versus scalar wave), but also dimensionality. We note that the acoustic scattering problem can be formulated in arbitrary 
spatial dimensions by replacing V <e R 3 with a bounded open set flgi , and accordingly substituting the integral kernel 



2 a -ik\r-r'\ 



k z e 



4n\r-r'\ 



with K (d \r,r'): 



kU +2)/2 H™_ m (k\r-r>\) 

4i(27r|r-r'|)W- 2 » /2 ' 



where H v denotes the vth order cylindrical Hankel function of the second kind. In particular, 4iK (2) (r,r') = k 2 H (^|r- 
r'|) =-2f^ 2 [ln(^|r-r'|)]/7T+const+0(|r-r'| 2 ln(^|r-r'|)) and 2iif (1) (r,r') = ^e"^ |r_r ' 1 . Abrute-force computation reveals 
that the integral kernels K (d \r,r') have finite Hilbert-Schmidt bounds for d - 1,2,3, so they induce the corresponding 
Hilbert-Schmidt operators ^ <d) for acoustic scattering problems. The electromagnetic wave scattering can be modeled 
in lower spatial dimensions d — 1,2 by modifying the Born equation into 



E r 



c (r) = ((I- X 'S id) )E)(r):=E(r)-x [ E(r')K w (r,r')d d r' - ^V V- f E(r')K (d \r,r')d 

Jn k z [ Jn 



V 



reOdl 



The last term in the above integral equation drops off for d - 1, so ^ (1) is a Hilbert-Schmidt operator, just as c r§' (1) . For 
d - 2, we may define 



(2), 



d ™(r,r>) = G™(r,r>) + Gf(r,r>) 



-In . 

7r k\r — r'\ 



and notice that the Hermitian operator induced by the integral kernel VV'g (2) (r,r') is already of Hilbert-Schmidt type. 
This is because 



(w'-VV 2) (r,r')=-(«'-V')ln 



k\r-r'\ 



Ton 



I _ri-(r-r') 
ds' 7r|r-r'| 2 



Ton 



ds' 



r,r' ed£2 



14 
is bounded, provided that the boundary curve dQ is smooth. Thus, in place of Eq. [7j we have 



0<<f If [( W -V) (2) (r,r')][(»'-V')fl (2) (r,r')]ds'ld 
J da (Jen I 



s < +00, 

leading to the conclusion that the quadratic polynomial ^ (2) (/ + 2^ (2) ) is a Hilbert-Schmidt operator. In summary, 
c|W)(/ + 2c|W)) rf -i j s a Hilbert-Schmidt polynomial of degree d for spatial dimensions d - 1,2,3. 

We have tabulated the different aspects of analytic behavior for three-dimensional electromagnetic (vector wave) 
and acoustic (scalar wave) scattering in Table HI and compare the dimension-dependence of the two types of scattering 
problems in Table HU It might be noted that despite the drastically different analytic properties between these two 
types of scattering problems, the scalar wave approximation is still a popular approach to the treatment of interactions 
between light and "objects much larger than wavelength", such as in the scalar diffraction theory [Ref. 6, pp. 478-482]. 
This is not particularly surprising as the overall spectral behavior of 'S can be very close to that of a Hilbert-Schmidt 
operator 7 in the short-wavelength limit (cf. the discussion of spectral series Y.\^ a t^) W 2 |l + 2A| 4 in 34.21 1. 

However, the goodness of approximation generated from scalar wave models of light-matter interaction cannot be 
taken for granted in all scenarios. We caution that for dimensions d — 2 and 3, there are some practically relevant 
regimes ("large wavelength" |20Ll2lll . "small particle" I22L 12311 ) where the salient contributions to the spectrum ct*(^) 
come from two sequences of points aggregating around the points {0} and {-1/2}. The latter sequence points to a 
universal resonance mode \l% - -1/2 (i.e e r - -1) that attracts scattering eigenvalues, which is a property intrinsic to 
electromagnetic scattering, unfound in scalar wave phenomena. 
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